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ABSTRACT 

Flapping-wing  propulsion  is  studied  experimentally  through  thrust  measurements 
and  flow  visualization.  The  objective  of  the  research  is  to  provide  further  insight  into  the 
aerodynamics  of  flapping-wing  micro  air  vehicles  (MAVs).  Experimental  work  is 
conducted  in  the  NPS  1.5  m  x  1.5  m  in-draft  wind  tunnel.  A  previously  constructed 
model  is  suspended  by  thin  wires  and  is  used  to  measure  the  thrust  performance  of  the 
flapping-wing  MAV  .  For  this  experiment,  the  model  is  tested  in  four  configurations; 
three  with  varying  wing  mount  stiffness  and  the  fourth  with  an  articulated  pitch 
mechanism.  Thrust  is  indirectly  determined  using  a  laser  range-finder  to  measure  stream- 
wise  displacement  of  the  model.  Three  methods  of  flow  visualization  are  attempted  to 
gain  further  insight  into  the  flow-field  around  the  MAV.  First  tufts  are  placed  on  and 
around  the  model  to  identify  the  flow-field.  Second,  a  smoke  rake  placed  outside  the 
tunnel  is  used  to  route  smoke  into  the  test  section.  Thirdly,  a  smoke  wire  system  is  used 
to  produce  smoke  in  the  test  section.  Experimental  results  are  compared  with  flow 
visualization  results  and  previous  experimental  and  numerical  work. 
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LIST  OF  SYMBOLS 


AR 

aspect  ratio,  b2/S 

b 

effective  wing  span 

b* 

total  wing  span 

c 

chord  length 

cd 

drag  coefficient  per  unit  span,  DI(qodc) 

cD 

drag  coefficient,  Dl(qc£)=  Cjb 

cP 

power  coefficient,  power/(g«»S Voo)=  -C,y-  Cma 
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thrust  coefficient,  TI(qo£) 

D 
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oscillation  frequency  in  Hz 

Jm 

pendulum  frequency  of  installed  swinging  model 

8 

acceleration  due  to  gravity 

ky 

vertical  plunge  amplitude  in  terms  of  c 

k 

reduced  frequency,  2nfclVoo 
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effective  pendulum  length 

q<=o 

free-stream  dynamic  pressure,  x/2  pJV^ 
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wing  area,  be 
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time 
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thrust 
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free-stream  velocity 

Xp 
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vertical  displacement  in  terms  of  c 

a 

angle  of  attack  (AOA) 

CCeff 

effective  angle  of  attack 

(Xgeom 

geometric  angle  of  attack 

&ind 

induced  angle  of  attack 
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sinusoidal  pitch  amplitude 

fa 

phase  difference  between  pitch  and  vertical  plunge 
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propulsive  efficiency,  CtICp 
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taper  ratio,  c/co 
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sweep  angle 

Poo 

free-stream  air  density 
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non-dimensional  time,  tVJc 
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circular  frequency,  27t/ 
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I.  INTRODUCTION 

A.  OVERVIEW 

In  this  paper  the  effects  of  wing  mount  stiffness  on  the  aerodynamic  performance 
of  a  flapping-wing  Micro  Air  Vehicle  (MAV)  was  investigated.  In  addition,  flow 
visibility  methods  were  utilized  in  an  attempt  to  gain  further  knowledge  of  the  flow  field 
about  the  MAV.  While  there  are  many  potential  propulsive  methods  for  MAVs, 
including  conventional  or  rotary  wing  designs,  this  paper  builds  upon  previous  research 
conducted  at  NPS  to  increase  the  knowledge  base  of  the  flapping-wing  propulsion  system 
for  the  MAV  application. 

The  experimentation  for  this  paper  was  conducted  in  the  NPS  1.5  m  x  1.5  m  in- 
draft wind  tunnel  using  a  previously  developed  parallel  flapping-wing  MAV.  To  obtain 
thrust  data,  the  MAV  was  suspended  with  wire  from  the  roof  of  the  tunnel  with  thrust 
indirectly  measured  using  a  laser  rangefinder  to  measure  the  streamwise  displacement  of 
the  MAV. 

B.  BACKGROUND 

With  a  shift  towards  more  diverse  military  operations  in  recent  years,  an 
operational  need  for  locally  owned  and  operated  reconnaissance  assets  has  been 
identified.  The  recent  development  of  new  technologies  such  as  micro-electromechanical 
systems  (MEMS),  tiny  CCD-array  cameras,  equally  small  infrared  sensors  and  chip-sized 
hazardous  substance  detectors  has  made  the  development  of  small  airborne 
reconnaissance  vehicles  viable.  To  this  end,  the  Defense  Advanced  Research  Projects 
Agency  (DARPA)  is  funding  research  into  the  development  of  MAVs.  A  MAV  is 
defined  by  DARPA  as  "an  affordable,  fully  functional,  militarily  capable,  flight  vehicle, 
limited  to  15  cm  in  length,  width  or  height."[Ref.  1].  These  vehicles  must  be  able  to  stay 
aloft  for  20  to  60  minutes,  cover  a  distance  of  up  to  10  km  and  carry  a  payload  weighing 
up  to  20  grams.  The  flight  regime  in  which  MAVs  will  fly,  is  for  the  most  part 


unstudied.  The  low  Reynolds  number  flight,  which  will  be  encountered  by  MAVs,  has 
been  up  until  this  point,  strictly  the  domain  of  insects  and  small  birds. 

Essential  to  the  success  of  a  MAV  design  is  the  requirement  for  high-density 
propulsion  and  power  sources  to  fulfill  the  MAVs  propulsive  requirements.  Previous 
research  indicates  that  a  parallel  flapping-wing  may  produce  higher  propulsive  efficiency 
than  a  conventional  propeller  driven  or  rotary  winged  vehicle  in  low  Reynolds  number 
flow.  The  research  summarized  in  this  paper  builds  upon  previous  research  conducted  in 
the  study  of  flapping-wing  MAVs. 


C.        FLAPPING- WING  PROPULSION 

In  1909  and  1912,  Knoller  [Ref.  2]  and  Betz  [Ref.  3]  respectively,  independently 
observed  that  a  flapping-wing  creates  an  effective  angle  of  attack,  ote,  resulting  in  a 
normal  force  vector  with  both  lift  and  thrust  components.  The  Knoller-Betz  effect  is 
illustrated  in  Figure  1. 
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Figure  1.  Effective  Angle  of  Attack  Produced  by  Flapping-Wing. 


In  1922,  Katzmayr  [Ref.  4]  demonstrated  the  Knoller-Betz  effect  experimentally 
by  measuring  an  average  thrust  with  a  stationary  wing  in  a  sinusoidally  oscillating  wind 
stream.  This  field  of  research  was  further  advanced  in  1924  when  Birnbaum  [Ref.  5  &  6] 
observed  the  conditions  that  lead  to  flutter  or  the  generation  of  thrust  and  suggested  a 
flapping-wing  as  an  alternative  to  the  conventional  propeller.  Over  the  next  few  decades, 
much  research  was  conducted  on  the  aerodynamics  of  flutter  and  gust  response. 
However,  it  wasn't  until  the  mid- 1930s  that  further  research  was  conducted  by  von 
Karman  and  Burgers  [Ref.  7],  to  explain  the  theory  behind  the  generation  of  thrust  or 
drag  in  a  flapping-wing.  For  a  more  detailed  history  of  flapping-wing  research  refer  to 
Refs.  8  &  10. 

More  recently,  numerical  and  experimental  work  was  conducted  by  Jones  and 
Platzer  [Ref.  8].  In  their  research,  they  found  that  the  two-airfoil,  opposed-plunge 
configuration  was  the  configuration  which  appeared  most  promising  for  the  MAV 
application  due  to  this  configuration's  relatively  high  propulsive  efficiency  and  the 
inherent  stability  resulting  from  the  opposed  plunge  motion.  In  this  application,  the 
lower  airfoil  acts  as  an  image  airfoil  in  a  ground-effect  analysis.  The  opposed  plunge 
configuration  is  shown  in  Figure  2. 
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figure  2.  Airfoil  in  Ground  Effect. 

When  operating  in  ground  effect,  downwash  velocities  are  decreased.    The 
ground  presents  a  boundary  condition  requiring  the  normal  component  of  velocity  to  go 
to  zero  altering  the  streamline  pattern  around  the  airfoil.  In  an  analytic  ground-effect 
analysis,  a  mirror  image  of  the  airfoil  is  placed  across  the  plane  of  symmetry.  The  vortex 
system  of  this  image  airfoil  is  opposite  the  original  system.  In  order  to  satisfy  the 
boundary  condition,  the  vertical  velocities  induced  by  these  two  image  vortex  systems 
will  cancel  at  the  plane  of  symmetry.  In  addition  to  vertical  velocity  going  to  zero,  the 
destructive  interference  of  the  counter-rotating  vortex  systems  diminishes  the  influence  of 
the  trailing  vortices  on  the  airfoil. 

The  parallel  flapping-wing  configuration  studied  in  this  paper  approximates  the 
beneficial  aspects  of  ground-effect  when  flying  out  of  ground-effect.  In  this  three- 
dimensional  case,  the  reduced  influence  of  the  tip  vortices,  lessen  the  induced  angle  of 
attack,  increasing  lift  and  thrust. 


II.    EXPERIMENTAL  ANALYSIS 


A.    EXPERIMENTAL  APPARATUS 


1. 


MAV  Construction 


The  MAVs  used  in  this  research  were  constructed  by  Dr  Kevin  Jones  [Ref.  10]. 
A  side  view  of  the  MAV  with  the  articulated  pitch  mechanism  can  be  seen  in  Figure  3. 
The  MAV  is  designed  around  an  exceptionally  small,  geared  stepping  motor  from  RMB 
Smoovy.  The  motors  are  5mm  in  diameter  with  the  25: 1  planetary  gear  system,  are 
approximately  25mm  long  and  weigh  a  mere  2.4g.  These  motors  produce  about  2.5  xlO"3 
Nm  of  torque  at  speeds  up  to  800  RPM. 


Figure  3.  MAV  with  Articulated  Pitch. 

The  brushless  motors  are  controlled  by  an  oscillatory  driver  circuit,  which  allows 

for  very  precise,  steady  and  reproducible  rotational  speeds.  As  the  control  circuit  and 

power  supply  are  large  and  bulky,  they  are  currently  mounted  external  to  the  model. 


Three-pole  power  is  fed  to  the  model  through  four  0.076mm  diameter  copper  wires  which 
support  the  model.  The  support  wires  connect  to  the  model  via  small  gold  plated  pins,  at 
the  nose  and  near  the  tail  of  the  model.  The  mass  of  the  wires  is  negligible  compared  to 
the  mass  of  the  model  and  they  are  flexible  enough  so  as  not  to  impede  the  model's 
motion. 

The  model's  body  has  composite  construction,  built  primarily  out  of  balsa  wood 
and  very  thin  graphite -epoxy  laminates.    The  motor  is  connected  by  thin-wall  silicon 
tubing  to  a  crankshaft  that  moves  the  flapping  beams  via  Scotch  yokes  constructed  of 
piano  wire.  The  wing's  leading  edge  is  a  teardrop  shaped  piece  of  balsa  wood  with 
carbon-fiber  ribs  running  to  the  trailing  edge  of  the  Japanese  tissue  wing  surface.  In  this 
paper,  the  wings  had  a  span  of  150mm  and  a  chord  of  35mm. 

The  wings  are  super-glued  to  the  flapping  beams  using  thin  carbon  fiber  strips 
whose  length  and  width  control  the  elasticity  of  the  joint.  The  elasticity  of  the  wing 
mount  allows  for  a  passive  feathering  mechanism.  The  wing  deflects  in  pitch 
proportionally  to  the  moment  about  its  leading  edge.  In  this  paper,  the  elasticity  of  the 
wing  joint  was  adjusted  by  decreasing  the  distance  between  the  wing's  leading  edge  and 
the  flapping  beam.  Three  different  distances  were  used  for  the  tested  model.  In  addition, 
an  active  pitch  control  was  added  to  the  model  by  gluing  arms  and  strings  to  the  model's 
body  and  wings  to  achieve  a  pure  plunge  motion  for  the  wing.  The  four  configurations 
were  tested  in  the  NPS  low-speed  in-draft  smoke  tunnel  and  the  results  compared. 

2.         Wind  Tunnel 

Experiments  were  performed  in  the  NPS  low-speed  in-draft  smoke  tunnel 
depicted  in  Figure  4.  The  tunnel  was  modeled  after  the  smoke  tunnel  of  the  Naval  Air 
Engineering  Lab  in  Philadelphia  [Ref.  12].  Air  is  ingested  from  inside  the  building 
through  a  square  4.5  x  4.5  m  inlet,  converging  through  a  9: 1  bell  shaped  contraction  to  a 
1.5  x  1.5  m  test  section.  Tunnel  speed  is  controlled  by  a  variable  pitch  fan  driven  by  a 
constant  speed  electric  motor.  Motor  and  fan  vibration  are  isolated  from  the  test  section 
by  rubber  sleeves  on  each  side  of  the  motor/fan  assembly.  The  tunnel  velocity  range  is  0 
to  9.5  m/s. 
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motor  with 
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Figure  4.  NPS  Low-Speed  In-draft  Wind  Tunnel. 

Late  in  1999,  Lt  Tim  Lund  rebuilt  the  neglected  wind  tunnel  [Ref.  11].  His 
extensive  repairs  included  fixing  the  fan,  resetting  the  variable  pitch  fan  blades  to  achieve 
lower  velocities,  removed  old  fixtures,  refinished  the  tunnel  floor  and  the  old  pitot-static 
system  was  replaced.  A  tunnel  calibration  was  performed  using  LDV  following  tunnel 
repairs  and  was  used  in  the  reduction  of  the  data  for  this  paper.  In  addition,  the  test 
section  turbulence  intensity  level  was  determined  and  remained  below  1.75%  with  tunnel 
velocities  over  1.5  m/s  [Ref.  11]. 

3.         Miscellaneous  Equipment 

A  pitot-static  tube  located  approximately  two  feet  below  and  one  foot  upstream  of 
the  model  was  used  to  determine  the  tunnel  velocity.  The  pitot-static  tube  was  connected 
to  a  MKS  Baratron  type  223B  differential  pressure  transducer.  The  transducer  voltage 
output,  which  is  linear  with  pressure,  was  read  by  a  Metex  3610D  digital  multimeter  that 
linked  the  voltage  data  to  a  PC  for  recording. 


A  frequency  strobe  light  was  used  to  set  the  flapping-wing  frequency.  The  strobe 
was  set  to  the  desired  frequency,  and  the  wing  frequency  adjusted  to  match.  The 
accuracy  of  the  strobe  was  checked  using  a  Monarch  Instrument  Tach  IV  Digital  Optical 
Tachometer.  Strobe  error  varied  with  frequency  to  a  maximum  of  1.5  percent  in  the 
frequencies  of  interest.  Actual  frequency  used  for  data  reduction  was  determined  after 
the  experiment  using  voltage  data  from  the  model  controller. 

A  Matsushita  Electric  Works  Model  ANL1651AC  laser  rangefinder  was  used  to 
determine  model  displacement  when  flapping.  The  sensor  accuracy  as  prescribed  by  the 
manufacturer  is  ±  lOOujn  ±  0.002  x  Ac  for  the  range  130mm  ±  35mm.  Its  measurable 
range  is  8  -  18  cm,  with  accuracy  decreasing  with  increased  distance  from  the  center 
point.  The  rangefinder  was  positioned  behind  the  model  at  its  center  point  distance  of  13 
cm.  During  runs,  this  corresponded  to  a  position  3  to  4  chord  lengths  behind  the  model 
and  was  not  thought  to  create  a  significant  flow  interference  effect.  The  output  of  the 
rangefinder  is  1  V/cm  displacement. 

A  Link  Instruments  DSO-2102  digital  storage  oscilloscope  was  used  to  measure 
voltage.  The  DSO-2102  was  connected  to  a  desktop  PC  for  display  and  recording  of  the 
voltage  signals.  The  oscilloscope  simultaneously  records  two  data  channels.  One 
channel  was  designated  for  the  laser  rangefinder,  and  this  data  was  subsequently  used  to 
determine  thrust.  The  second  channel  was  used  to  collect  controller  voltage  data  for 
eventual  determination  of  flapping  frequency  and  to  synchronize  data.  The  oscilloscope 
features  32  kilobyte  buffers  per  channel. 

A  Setra  EL410D  digital  balance  was  used  to  measure  the  mass  of  the  model 
before  each  run.  The  scale  has  an  accuracy  of  ±  0.00 lg.  The  model's  mass  was  found  to 
vary  by  as  much  as  0.005g  from  day  to  day,  due  to  humidity,  dust  or  other  environmental 
contamination.  A  Rosco  model  4500  smoke  generator  was  used  for  flow  visibility.  A 
smoke  rake  was  constructed  to  match  the  output  of  the  smoke  generator. 

A  Sony  DCR-VX1000  digital  Video  camera  was  used  to  capture  video  images  of 
the  MAV  in  motion.  Using  the  manual  focus,  auto  exposure,  a  shutter  speed  of  1/8  sec 
and  the  strobe  light  set  at  approximately  12  Hz  (just  slightly  faster  than  the  MAV 


flapping  frequency),  acceptable  video  images  of  the  MAV  were  obtained.  Capturing 
flow  visualization  images  using  the  smoke  rake  or  smoke  wire  proved  to  be  more 
difficult.  A  Kodak  DC  260  digital  camera  was  used  to  capture  flow  visualization  images 
(thread  tell-tales).  Best  results  were  achieved  by  setting  the  manual  focus  to  18  inches, 
the  exposure  to  0.5  sec  (the  slowest  manual  timed  exposure)  and  the  strobe  light  and 
MAV  flapping  frequency  set  to  approximately  8  Hz.  This  configuration  lead  to  some 
blurring  as  approximately  4  cycles  were  captured  per  image. 
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EXPERIMENTAL  CONFIGURATIONS 


Wing  Span 

b* 

150  mm 

Wing  Chord 

C 

35  mm 

Aspect  Ratio 

AR=b2/S=b/c 

4.29 

Taper  Ratio 

A=Ct/Co 

1 

Sweep  Angle 

A 

0° 

Velocity 

Vro 

0-6m/s 

Flapping  Frequency 

F 

8,  10,  12,  14  Hz  (approximate) 

Mean  Wing  Separation 

Aym 

0.8429c 

Plunge  Amplitude 

H 

25.4  mm 

Reduced  Frequency 

K 

0.29  -  °° 

Table  1.  General  Experimental  Configuration. 


Wing  Stiffness  (Low) 

K 

4.746  x  10"3  Nm 

Wing  Stiffness  (Medium) 

K 

6.384  x  10"3  Nm 

Wing  Stiffness  (High) 

K 

14.17  x  10"3Nm 

Wing  Stiffness  (Articulated  Pitch) 

K 

N/A 

Table  2.  Wing  Mount  Stiffness. 


C.         EXPERIMENTAL  METHOD 

Before  each  wind  tunnel  run,  the  model  with  mounting  pins  was  weighed  using 
the  Setra  scale.  The  model  was  then  hung  in  the  wind  tunnel  with  four  0.076mm 
diameter  copper  wires  and  mounting  pins  as  can  be  seen  in  Figure  5.  After  the  model 
was  installed  in  the  tunnel  the  laser  rangefinder  was  positioned  near  its  zero  voltage 
center  point,  13  cm  behind  the  model.  The  DSO  program  was  then  started  on  the  PC  and 
setup  to  record  data  at  1  KHz  (approximately  32  sec  of  data).  The  model  was  given  a 
push,  and  once  the  displacement  fell  within  the  limits  (500  mV/division,  10  divisions) 
displayed  in  the  DSO  program  (±  2.5  cm),  data  recording  began.  Three  separate  sets  of 
data  were  taken  before  either  the  tunnel  or  the  model  was  turned  on.  This  data  was  used 
to  determine  both  the  effective  pendulum  length  for  thrust  calculations  and  the  model's 
initial  position.  A  close-up  image  of  the  MAV  in  motion  with  laser  can  be  seen  in  Figure 
6. 


Figure  5.  MAV  Installed  in  Tunnel. 
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Figure  6.  MAV  in  Motion  and  Laser. 

Before  any  thrust  measurements  were  taken,  a  60-second  sample  of  pitot-static 
data  (sample  rate  of  1  Hz)  was  taken  to  establish  a  zero  velocity  voltage.  This  data  was 
recorded  on  the  PC  using  the  software  provided  with  the  Metex  3610D  DMM. 

Next,  a  set  of  thrust  measurements  were  taken  with  the  tunnel  off.  With  the 
controller  mode  selector  knob  in  the  variable  speed  position,  the  power  was  turned  on. 
The  default  power  setting  for  this  start-up  mode  was  set  so  that  the  model's  wings  would 
be  flapping  at  approximately  8  Hz.  The  strobe  was  set  to  8  Hz  to  verify  the  flapping 
frequency.  In  order  to  record  frequency  data,  the  TTL  output  from  the  controller  was 
connected  to  the  DSO.  With  the  25:1  gearbox  installed  in  the  model,  25  voltage  pulses 
from  the  controller  correspond  to  one  complete  cycle  for  the  model.  Once  the  model 
position  stabilized,  recording  of  the  position  data  by  the  DSO  program,  set  to  sample  at 
10  KHz,  was  started.  Once  again,  three  sets  of  data  were  recorded  at  each  frequency. 

The  strobe  frequency  was  then  increased  in  2  Hz  increments  up  to  14  Hz.  At  each 
2  Hz  increment,  the  signal  to  the  MAV  was  stepped  up  until  the  flapping  frequency 
matched  the  strobe  or  the  motor  stalled.  The  motor  tended  to  stall  if  the  load  became  too 
great  or  the  MAV  was  disturbed  by  an  external  disturbance. 
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Once  the  static  measurements  were  complete,  the  wind  tunnel  was  turned  on. 
Initially,  the  pitch  on  the  fan  blades  was  set  to  give  zero  tunnel  velocity.  The  pitch  was 
slowly  increased  incrementally  to  the  velocity  of  interest.  Once  the  velocity  stabilized, 
60-seconds  of  pitot-static  data  was  recorded  on  the  PC  (1  Hz  sample  rate). 

As  the  focus  of  this  paper  is  on  determining  the  feasibility  of  using  a  flapping- 
wing  mechanism  to  generate  thrust,  the  drag  of  the  MAV  is  required  so  that  the  gross 
thrust  generated  by  the  MAV  can  be  determined.  As  the  drag  of  the  MAV  varies  with 
wing  position  (maximum  drag  when  wings  are  together),  the  wings  were  flapped  at  an 
extremely  low  frequency,  on  the  order  of  0. 1  Hz,  to  determine  an  average  drag.  The  DSO 
was  used  to  record  the  position  data  from  the  laser  rangefinder  at  1  KHz  sample  rate. 
Again  three  sets  of  data  were  taken.  The  position  data  from  the  DSO  for  the  stiff  wing 
mount  run  at  approximately  5  m/s  and  0.1  Hz  can  be  seen  in  Figure  7. 

Stiff  Wing  Mount  Model  -  Raw  Drag  Data 

Flapping  Frequency  -  0.1 1  Hz 

Tunnel  Velocity  -  5.14  m/s 

250 


yV^     ^%^~/^    ^Vy**/"^ 


200 


5   150 

a. 


100 


50 


0         4000        8000        12000       16000       20000       24000       28000       32000 

Sample 

Figure  7.  Change  in  MAV  Drag  with  Wing  Position. 
Once  the  drag  data  had  been  recorded,  thrust  measurements  were  taken  as  per  the 
static  case  at  8,  10,  12  and  14  Hz  (if  the  motor  allowed).  Once  the  measurements  were 
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made,  the  tunnel  speed  was  increased  until  a  tunnel  speed  of  5  to  6  m/s  was  reached.  At 
tunnel  speeds  in  this  range,  the  model's  drag  displaced  it  sufficiently  to  force  the 
displacement  data  outside  the  display  limits  of  the  DSO  software.  With  tunnel  speeds 
from  0-5.9  m/s  and  flapping  frequencies  from  8  -  14  Hz  reduced  frequencies  in  the 
range  k  =  0.29  -  ~  were  observed. 

D.    EXPERIMENTAL  DATA  REDUCTION 

At  each  test  point,  three  different  sets  of  data  were  recorded;  voltage  from  the 
differential  pressure  transducer,  output  from  the  controller  and  voltage  data  from  the  laser 
rangefinder.  To  reduce  the  pressure  transducer  data  to  velocity  [m/s],  Equation  1, 
obtained  from  an  LDV-  calibration  completed  by  Lt  Lund  [Ref.  11],  was  used. 

WZ  =  0.4819(mV)0489  U) 

The  pressure  data  file  was  opened  in  Excel  and  converted  to  velocities  using 
Equation  1.  The  average  and  standard  deviation  of  the  60-second  sample  were  then 
taken.  To  account  for  changing  environmental  conditions,  an  initial  and  final  pressure 
transducer  reading  were  taken.  A  linear  change  in  the  zero  condition  was  assumed  and 
the  difference  between  the  initial  and  final  condition  was  divided  by  the  number  of 
velocities  and  added  to  each  velocity  reading.  The  test  velocities  ranged  from  0-5.9 
m/s.  To  determine  the  model's  effective  pendulum  length  for  thrust  calculations  and  the 
model's  initial  position,  the  32-second  (1  Hz  sample  rate)  static  swing  data  file  was 
opened  in  Excel  and  the  model's  displacement  was  plotted  versus  sample  number  (see 
Figure  8). 
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Figure  8.  MAV  Displacement  vs  Sample  #. 
Next,  the  data  was  edited  so  that  an  integer  number  of  cycles  were  represented. 
The  number  of  cycles  and  the  number  of  samples  within  the  cropped  data  set  were 
counted  and  the  frequency  of  the  hanging  model  was  determined.  The  installed  model's 
frequency,  fm,  was  approximately  1.07  Hz  for  the  configurations  studied.  The  model's 
pendulum  length,  L,  was  then  calculated  using  Equation  2  where  g  is  the  acceleration  due 
to  gravity. 

3—  (2) 


L  = 


(2xfJ 


The  model's  initial  position  was  determined  by  averaging  the  cropped  data  set. 
Next  a  MATLAB  program  (Appendix  A),  based  upon  one  written  by  Lt  Lund  [Ref.  1 1], 
was  written  to  reduce  the  raw  voltage  data  from  the  rangefinder  and  the  controller  and  to 
complete  error  calculations.  To  convert  the  voltage  data  to  model  displacement  in 
centimeters,  d,  Equation  3  was  used: 
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F  =  Mg  *  tan 


(  (d-d:\\ 

arcsm 


V 


v 


V*vpd*  10  24  (3) 

d  =  vpd*5.12-         P  K) 

256 

Where  vpd  is  the  volts  per  division  setting  in  the  DSO  software  display,  Vis  the 
voltage  output  from  the  laser  rangefinder,  dx  is  the  initial  displacement  of  the  model  in 
centimeters,  10.24  is  the  number  of  divisions  displayed  and  256  is  the  number  of  pixels 
displayed. 

The  force  required  to  displace  the  model  was  calculated  using  Equation  4: 

(4) 

L     )) 
Where  M  is  the  mass  of  the  model.  Equation  4  was  used  to  determine  both  the  MAV's 

drag  and  its  net  thrust  at  a  certain  tunnel  speed.  In  order  to  obtain  gross  thrust,  the  drag 

force  required  to  displace  the  MAV  model  was  subtracted  from  the  net  thrust. 

As  both  the  thrust  and  drag  of  the  MAV  vary  throughout  its  flapping  cycle,  the 
data  reduction  program  was  required  to  ensure  an  integer  number  of  cycles  were  used  for 
its  calculations.  To  accomplish  this,  the  voltage  data  from  the  controller  was  analyzed 
and  cropped  as  it  passed  through  a  certain  trigger  voltage.  One  voltage  cycle  is 
equivalent  to  a  motor  revolution  from  one  positive  peak  to  the  next  positive  peak.  One 
can  see  from  Figure  9  that  the  controller  gives  a  crisp,  repeatable  voltage  spike  for  each 
cycle.  As  the  gearbox  has  a  25:1  reduction  ratio,  25  peaks  equate  to  1  flapping  cycle  for 
the  MAV.  After  the  MATLAB  program  identifies  the  start  of  the  cycle,  it  counts  the 
number  of  times  the  controller  voltage  passes  through  the  trigger  voltage  (on  the  positive 
slope)  until  the  25th  time  when  the  program  calculates  the  voltage,  thrust,  and  frequency 
for  the  flapping  cycle. 

As  three  samples  of  data  were  taken  at  each  velocity/frequency  setting,  the  means 
from  each  sample  set  were  averaged.  The  reduced  data  was  saved  to  a  summary  file  for 
compilation  upon  completion  of  the  program.  The  program  also  allows  for  a 
comprehensive  summary  file  for  each  velocity/frequency  setting.  This  comprehensive 
summary  was  useful  for  trouble  shooting  during  construction  of  the  program. 
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Figure  9.  Controller  Voltage. 
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E. 


EXPERIMENTAL  RESULTS  AND  ANALYSIS 


1.         Experimental  Results 

In  this  paper,  the  effects  of  wing  mount  stiffness  on  the  performance  of  the 
flapping-wing  MAV  were  investigated.  Four  configurations  were  tested;  three  with 
increasing  wing  mount  stiffness  and  one  with  an  articulated  pitch  mechanism  which 
attempted  to  obtain  pure  airfoil  plunge. 

In  Figures  10  through  13,  the  thrust  obtained  from  each  configuration  is  plotted 
versus  velocity  at  approximately  8,  10,  12  and  14  Hz  flapping  frequency.  For  three- 
dimensional  plots  of  the  same  data,  see  Appendix  B. 

Low  Stiffness  Model 
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Figure  10.  Thrust  Performance  of  Low  Stiffness  Model. 
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Medium  Stiffness  Model 
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Figure  11.  Thrust  Performance  of  Medium  Stiffness  Model. 
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Figure  12.  Thrust  Performance  Fligh  Stiffness  Model. 
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Figure  13.  Thrust  Performance  of  Articulated  Pitch  Model. 

From  Figures  10  through  13  one  can  see  that  the  general  shape  of  the  thrust  vs 
velocity  plot  at  different  frequencies  is  basically  the  same  for  each  configuration. 
Initially,  the  model  produces  a  high  thrust,  which  quickly  reduces  to  a  minimum  as 
velocity  increases.  When  the  tunnel  speed  is  further  increased,  the  thrust  developed  by 
the  model  begins  to  slowly  recover.  One  can  see  that  as  the  flapping  frequency  increases, 
the  thrust- velocity  lines  move  upward  and  to  the  right  for  each  configuration. 

Figures  14  through  17  compare  the  thrust  performance  of  the  four  different 
configurations  at  each  frequency.  One  can  see  that  in  the  static  case  and  at  extremely  low 
tunnel  velocities,  the  models  with  the  more  flexible  wing  mounts  produce  greater  thrust. 
The  articulated  pitch  model  does  not  perform  well  in  the  static  and  low  speed  flow  cases. 
One  can  see  that  as  the  tunnel  velocity  increases  up  to  the  velocity  which  produces 
minimum  thrust,  the  thrust  performance  of  the  low  and  medium  stiffness  models  and  the 
high  stiffness  and  articulated  pitch  models  approach  each  other.  Once  past  the  velocity 
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that  produces  minimum  thrust  the  performance  of  the  medium  stiffness  and  articulated 
pitch  models  surpass  that  of  the  high  stiffness  and  low  stiffness  models. 
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Figure  14.  Configuration  Thrust  Comparison  @  8  Hz. 
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Frequency -10  Hz 
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Figure  15.  Configuration  Thrust  Comparison  @  10  Hz. 
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Figure  16.  Configuration  Thrust  Comparison  @  12  Hz. 
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Figure  17.  Configuration  Thrust  Comparison  @  14  Hz. 
From  the  above  results,  it  appears  that  the  best  wing  mount  configuration  for  a 
MAV  of  similar  construction  to  the  one  used  in  these  tests  would  be  one  that  had  a 
flexible  mount  at  low  velocities  and  a  high  stiffness  or  articulated  pitch  (pure  plunge) 
mount  at  higher  speeds. 

2.         Error  Analysis 

There  are  four  primary  sources  of  error  associated  with  the  experimental  thrust 
measurements.  These  four  include  the  error  due  to  the  average  deviation  of  thrust 
between  cycles,  laser  rangefinder  error,  oscilloscope  error,  and  weight/wire  length  error. 
The  sum  of  these  four  errors  is  presented  as  the  vertical  error  bars  in  Figures  10  through 
13.  As  the  in-draft  tunnel  vents  to  the  atmosphere,  the  velocity  measurements  taken  by 
the  pitot-static  system  were  highly  sensitive  to  wind  conditions.  Figures  10  through  13 
also  illustrate  the  fact  that  the  velocity  error  varies  greatly  day  to  day.    On  the  days  that 
the  data  for  the  low  stiffness  model  was  tested,  the  wind  was  relatively  calm.  However, 
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the  second  day  that  the  medium  stiffness  model  was  tested  was  extremely  windy  (> 
30mph).  The  resulting  error  in  velocity  was  as  much  as  22%.  The  velocity  error  bars 
displayed  in  Figures  10  through  13  represent  one  standard  deviation  added  and  subtracted 
from  the  mean  tunnel  velocity. 

The  first  source  of  thrust  error  is  thrust  average  deviation.  The  mean  thrust  for 
every  complete  flapping  cycle  of  the  model  was  calculated.  This  cycle  mean  thrust  of  the 
model  varied  from  cycle  to  cycle.  This  variation  was  due  to  a  variety  of  factors  but  the 
most  significant  contributors  were  that  the  tunnel  intake  was  indoors  and  that  the  outlet 
vents  to  the  atmosphere.  Although  the  tunnel  room  was  vented  to  the  outside  by  several 
banks  of  windows,  the  opening  and  closing  of  doors  and  other  windows  in  the  building 
caused  slight  velocity  drift  which  affected  model  displacement  and  thrust.  In  addition, 
variations  in  wind  velocity  also  affected  results.  Average  deviation  is  similar  to  standard 
deviation  in  that  it  provides  a  measure  of  the  width  of  the  data  about  the  mean.  It  is  a 
more  robust  estimator,  and  is  given  by  Equation  5  [Ref.  1]. 

1    N 

ADev(xl...xN)  =  —YJUj  -x]  (5) 

For  this  application,  x  represents  the  mean  of  each  cycle,  j  is  the  cycle  counter,  N  is  the 
total  number  of  cycles,  and  x  is  the  mean  of  all  cycles.  The  average  deviation  was 
calculated  at  every  frequency  and  velocity  and  is  considered  here  to  be  an  estimated  error 
bound.  It  is  the  predominant  source  of  error  for  velocities  above  approximately  4  m/s. 

The  second  source  of  thrust  error  was  the  laser  rangefinder.  Rangefinder  error 
consists  of  a  fixed  error  at  the  center  point  (13  cm).  A  linear  error  proportional  to  the 
distance  from  the  center  point  is  added  to  the  fixed  error  at  all  other  points.  For  the 
ranges  encountered  here,  the  rangefinder  accuracy  is  given  by  the  Equation  6: 

Range  Error —  0.  Imm  ±  0.2%  *  Distance  from  Center  Point  (6) 

The  third  source  of  thrust  error  was  due  to  the  digital  oscilloscope.  An  error  due  to 
the  maximum  resolution  of  the  oscilloscope  was  added  to  all  voltage  measurements.  This 
error  varies  with  the  selected  maximum  value  of  the  displayed  vertical  axis.  The  highest 
accuracy  is  obtained  when  the  voltage  signal  approaches  the  maximum  displayable,  but 
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does  not  exceed  the  display  window.  The  oscilloscope  accuracy  analysis  indicated  that 
the  actual  oscilloscope  error  was  always  less  than  half  of  its  resolution.  Statistically,  this 
is  reasonable  in  that  errors  would  tend  to  occur  on  both  sides  of  the  actual  voltage, 
offsetting  each  other.  In  the  interest  of  a  conservative  analysis,  however,  the  maximum 
resolution  was  used  as  the  error.  Due  to  the  requirement  to  capture  the  maximum 
displacement  of  the  model  at  tunnel  velocities  approaching  6  m/s,  the  DSO  resolution  had 
to  be  set  at  0.5  Volts  per  division.  Consequently,  the  error  associated  with  DSO 
resolution  was  the  greatest  contributor  to  thrust  error. 

The  fourth  source  of  error  was  that  associated  with  the  thrust  calculations.  As 
discussed  in  the  Experimental  Data  Reduction  section,  to  determine  the  thrust  produced 
by  the  model  the  weight  and  length  of  the  model's  wires  had  to  be  determined.  In  order 
to  account  for  the  uncertainties  in  these  measurements,  a  thrust  error  calculation  was 
carried  out  using  the  maximum  weight  and  minimum  wire  length  and  was  subtracted 
from  the  mean  thrust  for  each  data  set. 
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III.       FLOW  VISUALIZATION 

A.         FLOW  VISUALIZATION  METHODS 

Three  methods  of  flow  visibility  were  used  to  try  and  gain  insight  into  the  flow 
field  about  the  MAV.  Initially,  a  smoke  rake  and  fog  machine  located  outside  the  tunnel 
intake  were  used  to  feed  smoke  into  the  test  section.  Then  tell-tale  tufts  were  placed  on 
and  around  the  MAV's  wings.  Finally  a  smoke  wire  smoke  system  was  used  to  burn  fog 
fluid  in  the  test  section  near  the  MAV. 

1.         Smoke  Rake  and  Fog  Machine 

Initially,  attempts  to  document  the  flow  field  around  the  MAV  were  conducted 
with  a  smoke  rake  located  just  outside  the  tunnel  intake  connected  to  a  Rosco  4500 
smoke  generator  by  a  25  foot  long  six-inch  ID  flexible  duct.  A  smoke  rake  was 
constructed  with  an  eight  foot-long  main  tube  of  six-inch  ID,  schedule  40  PVC  pipe  with 
16-1.5  inch  ID  smoke  tubes  cemented  to  the  main  tube.  This  rake  was  mounted  on  a 
height  adjustable  stand  and  can  be  seen  in  Figure  18. 


Figure  18.  Smoke  Rake  In  Front  of  Tunnel  Inlet. 
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After  several  attempts  a  couple  of  problems  with  using  the  smoke  rake  located  in 
front  of  the  tunnel  intake  were  discovered.  As  the  fog  fluid  is  burned  to  generate  smoke, 
it  is  much  warmer  than  ambient  air.  Consequently,  when  the  smoke  enters  the  bottom  of 
the  rake,  it  tends  to  jet  past  the  first  few  smoke  tubes  and  most  of  the  smoke  exits  the  rake 
from  the  top  eight  tubes.  The  smoke  continued  to  rise  as  it  entered  the  tunnel's  9:1 
contraction. 

In  addition,  as  the  smoke  passed  through  the  contraction,  the  smoke  stream  tubes 
began  to  wrap  around  each  other,  forming  a  large  spiral  of  smoke  by  the  time  it  reached 
the  test  section.  As  the  wing  span  of  the  MAV  was  only  150mm  and  the  diameter  of  the 
smoke  spiral  was  approximately  30  cm  in  diameter  at  the  test  section,  the  smoke  totally 
enveloped  the  MAV  and  this  method  was  fairly  useless  for  flow  visibility  purposes. 
Consequently,  other  methods  of  flow  visibility  were  investigated. 

2.  Tell-Tale  Tufts 

As  Jones  [Ref.  10]  had  some  success  using  single  strands  of  cotton  thread  (a 
thread  is  usually  constructed  of  three  or  more  intertwined  strands)  glued  to  the  MAV's 
wings,  this  method  was  attempted.  To  reduce  the  strand  inertia  and  stiffness  problems 
associated  with  tufts  glued  to  the  wings,  single  strands  of  cotton  were  glued  to  a  balsa 
wood  spar  constructed  from  a  MAV  wing's  leading  edge.  Even  though  the  strands  were 
not  attached  to  the  rapidly  flapping-wing,  problems  with  the  inertia  and  stiffness  of  the 
fibers  were  encountered.  Once  an  optimal  fiber  length  was  determined  (approximately 
40  mm)  the  spar  was  placed  in  various  positions  around  the  MAV.  The  best  results  were 
obtained  with  the  spar  glued  to  the  center  rear  of  the  MAV's  fuselage  with  the  strands 
passing  between  the  two  flapping-wings. 

3.  Smoke  Wire  System 

The  third  method  used  to  document  the  flow  field  about  the  MAV  was  a  smoke 
wire  system.  For  this  method,  a  bare  strand  of  0.022"  nickel  chromium  (nichrome)  wire 
was  fixed  to  span  the  tunnel  just  ahead  of  the  MAV  in  the  test  section.  Electrical 
connectors  were  attached  to  the  end  of  the  wire  on  the  floor  of  the  tunnel,  and  on  an 
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adjustable  brass  tube  just  above  the  MAV.  The  adjustable  brass  tube  also  had  a  port  to 
which  a  clear  silicone  tube  was  attached  which  delivered  fog  fluid  from  a  pressurized  jug 
outside  the  test  section.  The  jug  was  pressurized  using  a  bike  pump  through  one  of  two 
adapters  bolted  through  the  jug's  cap.  The  clear  silicone  tube  was  attached  to  the  second 
adapter  and  fog  fluid  flow  was  controlled  using  a  pair  of  hemostat  clamps. 

Power  for  the  smoke  wire  was  provided  by  a  Variac  set  to  provide  15-20V.  The 
Variac  setting  was  dependent  upon  tunnel  speed  and  fog  fluid  flow  rate.  The  model  was 
installed  in  the  tunnel  on  a  sting  that  was  mounted  on  a  60  cm  span  symmetric  airfoil  that 
was  in  turn  mounted  to  a  two  axis  traverse.  The  flow  was  illuminated  using  a  Smith- 
Victor  Corp  750W  Photographic  Light. 

Different  configurations  of  the  wire  were  tried  to  determine  the  optimum  design. 
Initially,  loops  were  bent  in  the  wire  so  that  the  fog  fluid  would  pool  there  and  produce 
smoke  for  an  extended  period.  This  method  did  cause  the  fluid  to  pool  but  the  loops 
created  significant  vortex  shedding.  Next,  0.005"  bare  copper  wire  was  tied  around  the 
nichrome  wire  to  create  longer  duration  smoke.  Vortex  shedding  was  once  again 
observed,  however  the  magnitude  was  much  less  than  that  caused  by  the  looped  nichrome 
wire.  The  looped  wire  worked  fairly  well,  however  after  a  few  fog  cycles,  the  copper 
wire  knots  slid  down  the  nichrome  wire.  Finally,  a  straight  bare  nichrome  wire  was  used. 
Though  some  vortex  shedding  was  still  observed  at  higher  flow  speeds,  this  configuration 
gave  the  most  consistent  smoke  plane.  The  smoke  wire  set-up  can  be  seen  in  Figures  19 
and  20. 
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figure  ZU.  Smoke  wire  set-up  in  lunnei. 
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Flow  images  from  the  smoke  wire  technique  were  captured  using  a  Sony  DVC- 
2000  digital  video  camera  and  a  Kodak  DC  260  digital  camera.  When  using  the  Kodak 
camera,  the  focus  was  manually  set  at  either  24  or  36  inches  and  the  flash  was  turned  on. 


B. 


FLOW  VISIBILITY  RESULTS 


1. 


Tell  -Tale  Tufts 


Using  the  tuft  method,  the  best  results  were  obtained  when  the  tuft  spar  was 
placed  just  ahead  of  the  leading-edge  plane  of  the  wings  mounted  on  the  MAV  fuselage. 
In  Figure  21,  the  medium-stiff  winged  MAV  was  flapping  at  12  Hz  with  no  wind.  The 
MAV  was  lit  by  a  strobe  light  synchronized  with  the  model  and  the  image  captured  using 
the  Kodak  DC  260  digital  camera  with  the  focus  set  manually  to  18  inches  and  the 
exposure  set  to  0.5  sec. 


Figure  21.  Tell-Tale  Tufts,  Medium-Stiff  Wing,  12Hz  Static. 

With  the  given  frequency  and  exposure  time,  approximately  six  cycles  were 
captured  in  this  image.  This  can  clearly  be  seen  in  the  multiple  images  of  the  fourth 
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strand  from  the  right.  These  multiple  images  of  this  strand  indicate  that  the  flow  between 
the  wings,  near  the  centerline  of  the  MAV  is  quite  unsteady.  This  image  shows  the  three 
dimensionality  of  the  flow  field  around  the  MAV  as  the  flow  is  entrained  as  the  wings 
move  apart.  Figure  22  illustrates  the  outward  flow  from  between  the  wings  as  they  move 
together.  This  image  is  of  the  MAV  with  the  medium-stiff  wing  mount,  flapping  at  12 
Hz  with  flow  speed  of  approximately  2.9  m/s. 


Figure  22.  Tell-Tale  Tufts,  Medium-Stiff  Wing,  12Hz,  2.9  m/s  Wind  Speed. 

2.         Smoke  Wire  System 

Initially,  attempts  were  made  to  illuminate  the  smoke  flow  using  the  strobe  light. 
It  quickly  became  clear  that  the  output  of  the  strobe  would  not  be  sufficient  for  capturing 
smoke  images.  In  addition,  as  there  were  no  means  available  to  synchronize  the  camera 
with  the  strobe,  an  exposure  length  longer  than  one  complete  cycle  was  necessary  to 
ensure  that  the  strobe  flash  was  captured.  As  the  flow  around  the  MAV  is  unsteady,  these 
attempts  to  capture  an  image  over  more  than  one  cycle  lead  to  significant  blurring  of  the 
smoke  plane.  Finally,  a  Smith- Victor  Photographic  Light  was  used  to  illuminate  the  flow 


30 


and  the  Kodak  DC260's  flash  was  turned  on.  The  images  in  Figures  23  through  26  were 
taken  from  approximately  three  to  four  feet  downstream  of  the  MAV  with  the  manual 
focus  set  to  two  or  three  feet. 


Figure  23.  MAV  Wing  Tip  Vortex  Induced  Wake  Structure. 

In  Figure  23,  the  smoke  wire  is  positioned  approximately  six  inches  upstream  of 
the  MAV's  port  leading  edge.  One  can  clearly  see  the  wake  structure  developed  by  the 
wing  tip  vortices  as  they  are  shed  during  the  flapping  motion.  Looking  upstream,  as  the 
upper  and  lower  wings  move  apart,  they  respectively  shed  counter-clockwise  and 
clockwise  rotating  vortices.  Correspondingly,  when  the  upper  and  lower  wings  move 
together,  they  respectively  shed  clockwise  and  counter-clockwise  rotating  vortices.  The 
wing  tip  vortices  can  be  more  clearly  observed  in  Figure  24. 
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Figure  24.  Close  Up  of  MAV  Wing  Tip  Vortices. 

The  vortices  that  are  clearly  visible  in  Figure  24,  were  formed  as  the  upper  and 
lower  wings  moved  together.  As  these  vortices  grew  and  coalesced,  they  formed  a  jet, 
which  flowed  outwards  away  from  the  MAV's  centerline.  The  interaction  of  these 
vortices  also  pushed  the  tip  vortices  outwards.  Once  the  wings  reached  their  minimum 
separation,  they  stopped  shedding  their  vortices  and  the  smoke  plane  remained  flat  as  can 
be  seen  in  the  vertical  smoke  line  just  outboard  of  the  wing  tips  joining  the  two  vortices. 
When  this  image  was  taken,  the  MAV's  wings  had  just  begun  to  move  apart.  The  tip 
vortices  which  were  just  forming,  are  not  visible  but  the  effects  of  the  inward  flowing  jet 
they  formed,  can  be  seen  in  the  deflected  smoke  plane  between  the  wings  to  the  right  of 
the  smoke  wire  in  the  image. 

The  alternating  vortices  of  the  upper  wing  and  their  effects  on  the  smoke  plane 
can  be  seen  in  Figure  25  taken  from  above  and  downstream  of  the  MAV.  In  this  image 
the  clockwise  rotating  vortex  nearest  the  MAV  was  shed  from  the  upper  wing  as  the 
wings  moved  together.  The  effects  of  the  outward  flowing  jet  created  by  this  vortex  set 
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on  the  upper  portion  of  the  smoke  plain  can  be  seen  in  the  deflection  of  the  top  of  the 
plain  to  the  left  (seen  to  the  left  of  the  vortex  in  this  image). 


Figure  25.  Alternating  Vortices  Shed  from  the  MAV's  Upper  Wing. 

The  second  vortex  aft  of  the  MAV  is  counter-clockwise  rotating  and  was  shed 
from  the  upper  wing  as  the  wings  moved  apart.  The  effects  of  the  inward  flowing  jet 
created  by  this  second  vortex  set  on  the  upper  portion  of  the  smoke  plain  can  be  seen  in 
the  deflection  of  the  top  of  the  plain  to  the  right  (seen  above  and  to  the  left  of  this  vortex 
in  the  image). 
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The  final  image  obtained  by  the  smoke  wire  method  can  be  seen  in  Figure  26. 
This  image  was  taken  with  the  smoke  wire  positioned  in  front  of  the  MAV's  nose  directly 
along  its  centerline.  In  this  image  one  can  see  the  effects  of  the  boundary  layer  over  the 
MAV's  nose  in  the  upward  deflection  of  the  smoke  plane.  As  this  image  was  taken,  the 
wings  were  at  their  minimum  separation;  consequently,  the  motion  of  the  wings  is 
approaching  zero  and  the  smoke  plane  is  relatively  undisturbed  above  the  upper  wing. 
Just  behind  the  upper  wing,  the  smoke  plane  is  displaced  upwards  by  the  clockwise 
rotating  (as  seen  from  this  side  of  the  smoke  plane)  vortex  shed  from  the  upper  wing  as 
the  wings  moved  together.  As  expected  in  Figure  26,  the  flow  is  relatively  two 
dimensional  when  compared  to  the  flow  about  the  MAV's  wing  tips. 


Figure  26.  Two  Dimensional  Flow  Over  the  MAV's  Centerline. 

Figures  27  through  35  summarize  of  the  flow  field  information  obtained  using  the 
smoke  wire  system.  In  Figure  27,  the  MAV's  wings  are  at  their  minimum  separation  and 
the  vortices  shed  during  the  closing  of  the  wings  have  been  convected  downstream. 
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Figure  27.  MAV  Flow  Field  -  Wings  at  Min.  Separation  -  Steady. 
In  Figure  28,  the  MAV's  wings  are  beginning  to  move  apart.  The  wings  are 
deflected  due  to  aerodynamic  loading  and  the  wings'  inertia.  Due  to  the  movement  of  the 
wings  and  the  resulting  change  in  lift  generated  by  this  movement,  vortices  are  formed 
off  the  wings'  trailing  edges  and  wing  tips. 
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^igure  28.  MAV  Flow  Field  -  Wings  at  Min.  Separation  -  Moving. 
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In  Figure  29  the  wings  have  moved  apart  and  accelerated  to  their  maximum 
plunge  speed  at  their  median  separation.  As  the  wing  tip  vortices  grow  they  begin  to 
coalesce  and  form  a  jet  between  them,  flowing  towards  the  MAV's  centerline.  The 
interaction  of  the  vortices  also  pushes  the  vortex  pair  towards  the  MAV's  centerline.  The 
wind  convects  the  vortices  downstream. 
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Figure  29.  MAV  Flow  Field  -  Wings  at  Median  Separation. 
In  Figure  30  the  wings  are  decelerating  as  they  reach  their  maximum  separation. 
The  intensity  of  the  shedding  votices  decreases  and  they  continue  to  be  convected 
downstream. 
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"igure  30.  MAV  Flow  Field  -  Wings  Approaching  Max.  Separation. 
In  Figure  31  the  MAVs  wings  are  at  their  maximum  separation  and  are  no  longer 
moving.  Consequently,  the  wings  cease  to  shed  their  vortices. 
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^igure  31.  MAV  Flow  Field  -  Wings  at  Max.  Separation  -  Steady. 
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In  Figure  32  the  MAV's  wings  are  beginning  to  move  together.  One  can  see  that 
the  wings  deflect  in  the  opposite  direction  and  the  vortices  are  rotating  in  the  opposite 
sense  to  those  in  Figure  28. 
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^igure  32.  MAV  Flow  Field  -  Wings  at  Max.  Separation  -  Moving. 
In  Figure  33  the  wings  have  moved  apart  and  accelerated  to  their  maximum 
plunge  speed  at  their  median  separation.  As  the  wing  tip  vortices  grow  they  begin  to 
coalesce  and  form  a  jet  between  them,  flowing  away  from  the  MAV's  centerline.  The 
interaction  of  the  vortices  also  pushes  the  vortex  pair  away  from  the  MAV's  centerline. 
The  wind  convects  the  vortices  downstream. 
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"igure  33.  MAV  Flow  Field  -  Wings  at  Median  Separation. 
In  Figure  34  the  wings  are  decelerating  as  they  reach  their  minimum  separation. 
The  intensity  of  the  shedding  votices  decreases  and  they  continue  to  be  convected 
downstream. 
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'"igure  34.  MAV  Flow  Field  -  Wings  Approaching  Min.  Separation. 
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Finally,  in  Figure  35  the  MAVs  wings  are  at  their  minimum  separation  and  are  no 
longer  moving.  Consequently,  the  wings  cease  to  shed  their  vortices,  which  have  been 
convected  downstream. 
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^igure  35.  MAV  Flow  Field  -  Wings  at  Min.  Separation  -  Steady. 
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IV.       CONCLUSIONS 

With  a  shift  towards  more  diverse  military  operations  in  recent  years,  an 
operational  need  for  locally  owned  and  operated  reconnaissance  assets  has  been 
identified.  The  recent  development  of  new  technologies  has  made  the  construction  of 
small  airborne  reconnaissance  vehicles  viable.  Nature  hints  that  flapping  wings  are 
efficient  at  the  low  Reynolds  numbers  in  which  these  MAVs  would  operate. 
Consequently,  the  potential  efficiency  of  flapping  wings  as  a  means  of  MAV  propulsion 
makes  investigation  worthwhile. 

Experimental  testing  of  a  previously  constructed  MAV  was  conducted  in  the  NPS 
1.5  m  x  1.5  m  in-draft  wind  tunnel.    For  this  paper,  four  different  wing  configurations 
were  investigated;  three  with  varying  degrees  of  wing  mount  stiffness  and  one  with  an 
articulated  pitch  mechanism  which  attempted  to  obtain  a  pure-plunge  motion  for  the 
wings.  Thrust  was  indirectly  determined  by  measuring  the  stream-wise  displacement  of 
the  suspended  model  due  to  flapping.  The  steady-state  drag  was  subtracted  from  the 
computed  thrust  in  order  to  isolate  the  thrust  performance  of  the  MAV  and  minimize 
viscous  boundary-layer  drag  effects. 

After  collection,  reduction  and  analysis  of  experimental  data,  it  became  apparent 
that  in  order  to  optimize  the  thrust  performance  of  a  MAV  some  sort  of  active  pitch 
control  would  be  beneficial.  At  low  flow  speeds,  a  flexible  wing  would  produce  more 
thrust,  while  at  higher  speeds  a  stiffer  mount  or  active  pitch  control  to  approach  a  pure 
plunge  motion  would  be  more  effective. 

Images  collected  using  the  smoke  wire  method  were  extremely  beneficial  in 
determining  characteristics  of  the  flow  field  about  the  MAV.  Further  flow  visibility 
efforts  need  to  be  made  at  different  flapping  frequencies  and  tunnel  velocities  so  that  the 
flow  mechanics  responsible  for  the  changing  thrust  performance  can  be  better 
understood. 
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V.  RECOMMENDATIONS 

In  order  to  move  forward  in  the  development  of  the  flapping-wing  as  a  propulsion 
system  for  MAVs,  more  effort  must  be  put  into  visualization  of  the  flow  field  about  the 
MAV.  Initial  attempts  using  the  smoke  wire  system  were  fairly  effective  at  capturing 
images,  but  the  method  can  be  improved.  First,  a  constant  temperature  circuit  employed 
in  the  smoke  wire  system  would  simplify  the  smoke  generating  process.  Secondly,  the 
use  of  a  laser  sheet  would  allow  images  of  the  cross  section  of  the  wake  behind  the  MAV 
to  be  captured.  Thirdly,  a  higher  intensity  strobe  light  or  multiple  synchronized  strobes 
may  allow  photographic  equipment  to  capture  "still"  images  of  the  flowfield.  This  would 
allow  the  photographer  more  control  over  which  portion  of  the  cycle  is  captured  on  film. 
Finally,  a  flow  highlighting  system  originating  from  the  MAV  wing  would  allow  for 
further  investigation  into  the  dynamic  stall  phenomenon  occurring  on  the  MAV  wings 
throughout  the  flapping  process  and  also  highlight  the  formation  of  the  wing  tip  vortices. 
Perhaps  the  use  of  titanium  tetra-chloride  in  the  in-draft  tunnel  could  be  investigated. 
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APPENDIX  A.  MATLAB  PROGRAM  FOR  DSO-2102  OSCILLOSCOPE  DATA 

REDUCTION 

iCe.pt    Sean  Duggan  -  MAV  •'  DSO  Data  Reduction  Program 

o 

%Data  Analysis  Process: 

%-  Weigh  model  with  pins  before  each  set  of  runs  (enter  weight  &  max  weight  [N]  below) 

%-  Install  model  in  tunnel,  center  laser,  swing  model  to  get  period  data  (1  KHz  sample 

rate) 

%-  Reduce  period  data  to  get  wire  length  [cm],  enter  below  (L= (g) / (2*pi*f req) ~2) 

%-  Use  period  data  to  determine  initial  model  displacement  [pixels] ,  enter  below 

%  Record  initial  velocity  (pressure)  data  before  tunnel  on 

%  Record  velocity  at  beginning  of  each  run 

%-  Record  final  velocity  with  tunnel  off 

%NOTE  -  True  velocity=0  -  4319*  (millivolts)  NC  .  489 

%        (From  LDV  /  Pitot  systerr  calibration  from  Lt  Tim  Lund's  Thesis) 

% 

%IMPORTANT:  Save  data  in  following  format: 

%       -  freqilOO)  ,  vel(vOOO)  ,  data  setd  thru  3) 

%         (ie.  lOOvOOQl,  100v0002,  lC0v00C3,  lOlvOOQl, , lxxv0003  etc.  for  each  vel) 

%       -  use  freq  100  for  drag  data  (IKHz  sample  rate) 

%       -  use  freq  101... lxx  for  thrust  data  (lOKHz  sample  rate) 

%       -  take  3  sets  of  data  at  each  freq/vel  setting 

%       -  save  data  at  each  velocity  in  its  own  directory  (ie.  vOOO,  vOOl,  v0C2...) 

& 

o 

%-  Archive  data  on  CD-R,  go  to  computer  with  MATLAB  to  reduce  data 

%-  After  data  is  reduced  at  each  velocity,  copy  data  in  Kean_Thrust.wri  file  to  summary 

file 

9- 

%KOTES:    If  you  open  Mear._Thrust.wri  in  EXCEL,  copy  data  to  summary  file  then  close 

%         Mean_Thrust.wri.   If  you  don't  close  the  file,  this  program  will  crash  when 

%         it  tries  to  delete  Mean_Thrust.wri  file 

% 

%-  When  data  reduction  complete,  subtract  drag  (IGOvXXXX)  from  thrust  at  each  vel  to  get 

%   corrected  thrust 

%-  Compile  velocity  data  in  one  summary  file 

%-  Determine  avg  vel  and  standard  deviation  of  vel  and  copy  to  thrust  summary  file 

clear 

clear  all 

delete  H: \docs\Results\Mean_Thrust.wri; 

?.  ************** ******  R"^ii  ir}  this  *"*  a  *"~  t  ****************************** 

filelocation='F: \' ;       %  CHANGE  FOR  EACH  COMPUTER  -  CD  Rom  drive  location 

date='17  May  00\';     %  CHANGE  FOR  EACH  RUN      -  Date  directory 

dirname='v01S\ ' ;  %  CHANGE  FOR  EACH  VELOCITY   Velocity  Directory 

filename= 'v018' 

filetype=  '  . prr. ' 

startf ilenumber=100;      %  First  file  number  -  Drag  data  (0.1Hz  flap) 

lastfilenumber=104;       %  CHANGE  FOR  EACH  VELOCITY  -  Last  file  number 

vol tsperdiv=0 . 5 ; 

weight=0 . 21329;  %  Model  weight  in  Newtons 

weight_max=0 .21334;  %  Max  weight  for  error  calcs 

wire_length=21. 8087386;  %  Cable  length  in  cm 

wire_length_min=21 . 8090208;    %  Min  wire  length  for  error  calcs 

disp_intial_pix=123 .0946956;   %  Initial  displacement  at  vO  fO  in  pixels 

disp_initial= (voltsperdiv*5 . 12) - (disp_intial_pix) * (voltsperdiv*10 .24/256) ; 

fidl  =  fopen{ 'H:\docs\Results\Mean_Thrust.wri' , 'w' ) ; 

fprintf (fidl, 'Filename\tVel\tVel\tMean\tCorrected\tTotal  Error \t ' ) ; 
fprintf (fidl, 'Avg  Dev\tLaser\tDSO\tWt  Wire\tMean\tAvg  Dev\tMin\tMax\n' ) ; 
fprintf (fidl, '\t\tStd  Dev\tThrust\tThrust\t+/  \tError\tError\tError\t ' ) ; 
fprintf (fidl, ' Error \t?requency\tFrequency'\tFrequency\tFrequency\n' ) ; 
fprintf  (fidl,  '\t[m/s]\t[m/s]  \t[N]\t[N]  \t  [N]  \t[N]\t[N]\t{N]\t[N]\t' ) ; 
fprintf (fidl, ' [Hz] \t [Hz] \t [Hz] \t [Hz] \n\n' ) ; 
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fclose(fidl) ; 


i=startf ilenumber : lastf ilenumber 
:or  k=l:3 

clear  data  %  Clears  dat 

if  i==100 

sample_rate=1000;  %  Sample 
else  sample_rate=10000;  %  Sample 
end 


matrix  eacr.  iteration 


rate 
rate 


in  riz  tin 

in  Hz  fm 


DSO 
DSO 


(0.1Hz  flap 
(Thrust) 


Drag; 


%  **********************  L0ads  data  *****************-********-***** 
disp( [' *******  Computing  ' ,num2str (i) , filename, num2str (k) , filetype, '  **' 
rawdata=load  ( [filelocat ion, date, dirname,num2str (i) , filename, nura2str(k) 
if  k  ==1 

datal ( : , 1) =rawdata ( : , 2 ) ;  %  CH  3  is  loaded  to  datal  col  1 

dataK : , 8)=rawdata( : , 1) ;  %  CH  A  is  loaded  to  datal  col  S 

%  ********  Converts  data  to  voltage  and  distance  (lV=lcm)  **i*****4-* 

dataK  :  ,  2)  =  (voltsperdiv*5.12)  -  (  (dataK:  ,1)  )*  (voltsperdiv*10. 24/256)  ) 

dataK : , 9 ) = (voltsperdiv*5 . 12 ) - ( (dataK : ,8) )* (vol tsperdiv*10 . 24/256) ) 

%  ********  Does  Nothing  ************************************** 

datal ( : , 3 ) =datal ( : , 2 ) ;  %  Does  nothing 

%  ********  converts  voltage  to  Force  [N]  ************************** 

datal ( : , 5) =  weight *  tan (asin ( (datal ( : , 3) -disp_initial) /wire_length) ) ; 


****']) 
filetype] 


%  ********  Computes  the  avg  dev  of  thrus 
%  DSO  data  variables 
if  i==100 

sig_trip  =  132; 
else  sig_trip  =  140; 
end 

sig_oldl  =  257; 
sig_old2  =  257; 
sig_old3  =  257; 
%  Voltage  variables 
cycle_volts=0; 
cycle_mean_volts=0 ; 
vo 1 1 s_suml  =  0 ; 
%  Thrust  variables 
cycle_thrust=0; 
cycle_mean_thrust=0 ; 
thrus  t_suml  =  0  ,- 
%  Frequency  variables 
eye 1 e_mean_f  req=  0 ; 
f req_suml=0; 
freql=0; 
x=l; 

%  Counter  variables 
nr_cyclesl=-l; 
cycle_data_pts=0; 
total_data_ptsl  =  0  ,- 
y=-l; 

data_pointsl=length (datal ( : , 1) ) 
for  j=l :data_pointsl   %  Looks  a 


freq  in  data 


Kin  data  value 

from  controller  (0.1  Hz  flap  - 

Min  data  value 

from  controller  (Thrust  runs) 

sig_oidl,  2,  a: 

id  3  are  the  3  data  points  prio 

to  the  current 

pt.  Use  257  if  the  trip  signal 

is  neg  voltage 

Use  -1  if  trip  is  a  +  voltage 

Drag  run) 


%  Sum  of  voltage  for  one  cycle 
%  Mean  volts  of  one  cycle 
>  Sum  of  voltages  of  all  cycles 

%  Sum  of  thrust  for  one  cycle 
%  Mean  thrust  of  one  cycle 
%  Sum  of  mean  thrust  of  all  cycles 

Mean  frequency  of  one  cycle 

Sum  of  frequencies  of  all  cycles 

%  Variable  for  min  frequency 

%  Kin  freq  counter 

%  Counter  for  total  number  of  cycles  in  data 
Number  of  data  pts  within  one  cycle 
Number  of  data  pts  of  all  cycles 

%  Counter  to  determine  cycles  (25  trips 

%  Counts  rows  of  data  matrix 

all  data  points 


1  cycle; 


(+/-  5V)  **' ) 


if  abs (datal (j ,2) )  >  5 

error  ('**  Voltage  is  outside  rangefinder 
end 
if  nr_cyclesl  >=  0  %  Ignores  partial  data  before  1st  cycle 

cycle_volts  =  cycle_volts  +  datal ( j -1, 2) ;  %  Sums  cycle  volts 

cycle_thrust  =  cycle_thrust  +  datal ( j-1, 5) ;  %  Sums  cycle  thrust 

cycle_data_pts=cycle_data_pts+l;  %  Counts  data  pts  in  cycle 
end 

%  Is  microswitch  tripped? 

if  datal ( j , 8) >=sig_trip  & 

sig_old2<sig_trip  & 

%  trip  voltage.  Change 

y=y+l; 

if  nr_cyclesl  ==  -1 
nr_cyclesl  =  0; 

end 


sig_oldl<sig_trip  &  ... 
sig_old3<sig_trip  %  This  assumes  e 
>,  <  signs  if  trip  is  a  +  voltage. 


neg 
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Pulses-  cycle 
Resets  count 


if  y  ==  25 

nr_cyclesl=nr_cyclesl+l;    % 

y=0; 

%  Volts 

cycle_mean_volts=cycle_volts/cycle_data_pts;  %  Cycle  mean 

datal ( j -1, 4) =cycle_mean_volts;  %  Stores  cycle  mean  volts 

volts_suml=volts_suml+cycle_mean_volts;  %  Sums  total,  volts 

%  Thrust 

cycle_mean_thrust=cycle_thrust/cycle_data_pts;  %  Cycle  sear. 

datal ( j -1 , 6) =cycle_mean_thrust;  %  Stores  cycle  mean  thrust 

thrust_suml=thrust_suml+cycle_mean_thrust;  %  Sums  tot  thrust 

%  Frequency 

cycle_mean_f req=sample_rate/cycle_data_pts;  %  Cycle  mean 

datal ( j-1, 7) =cycle_mean_freq;  %  Stores  cycle  mean  frequency 

f req_suml=freq_suml+cycle_mean_f req;  %  Sums  all  cycle  freqs 

if  datal (j-1, 7) >0 

freql(x)  =  cycle_mean_f req;   %  Save  nonzero  freqs 

x=x+l; 
end 

total_data_ptsl=total_data_ptsl+cycle_data_pts; 
%  Reset  variables 

%  Reset  cycle  voltage  for  new  cycle 
Reset  thrust  summation  for  new  cycle 
%  Reset  free  summation  for  new  cycle 
Reset  nr  of  data  ots  for  nev;  evele 


cycle_volts=0; 
cycle_thrust  =  0  ; 
cycle_f  req=0  ,- 
cycle_data_pts=0 ; 
end 
end 

sig_oldl  =  datal (j, 8) ; 
if  j>=3 

sig_old2  =  datal (j-1 
sig_old3  =  datal (j -2 
end 


Last  signal 


ilue 


%  The  data  value  before  the  last  one 
%  Three  data  points  back 


if  k  ==2 

data2 ( : , 1) =rawdata ( : , 2) ;  %  CH  3  is  loaded  to  datal  col  1 

data2 ( : , 8) =rawdata ( : , 1) ;  %  CH  A  is  loaded  to  datal  col  8 

%  ********  converts  data  to  voltage  and  distance  (lV=lcmj  ********** 

data2(: ,2)=( (voltsperdiv*5 . 12) - ( (data2( :,1) )* ( vol tsperdiv*10 .24/256) ) ) ; 

data2( : ,9)= ( (voltsperdiv*5 . 12) - ( (data2( :,8) )* ( vol tsperdiv* 10 .24/2  56)  )  )  ; 

%  ********  Does  Nothing  ************************************** 

data2 ( : ,3) =data2 ( : ,2) ;  %  Does  nothing 

%  ********  converts  voltage  to  thrust  [N]  ************************** 

data2 ( : , 5) =  weight * tan (asin ( (data2 ( : , 3) -disp_initial) /wire_length) ) ; 

%  ********  Computes  the  avg  dev  of  thrust  &  freq  in  data  *********** 

%  DSO  data  variables 


if  i==100 

sig_trip  =  132; 
else  sig_trip  =  140; 
end 

sig_oldl  =  257 
sig_old2  =  257 
sig_old3  =  257 
%  Voltage  variables 
cycle_volts=0; 
cycle_mean_volts=0 ; 
volts_sum2=0; 
%  Thrust  variables 
cycle_thrust=0 ; 
cycle_mean_thrust=0 ; 
thrus  t_sum2  =  0 ; 
%  Frequency  variables 
cycle_mean_freq=0; 
f req_sum2=0; 
freq2=0; 
x=l; 


%  Kin  data  value  from  controller  (0.1  Hz  flap  -  Drag  run) 
%  Kin  data  value  from  controller  (Thrust  runs) 

%  sig_oldi,  2,  and  3  are  the  3  data  points  prior 
%  to  the  current  pt .  Use  2  57  if  the  trip  signal 
%  is  neg  voltage.  Use  -1  if  trip  is  a  +  voltage. 


%  Sum  of  voltage  for  one  cycle 
%  Mean  volts  of  one  cycle 
'<;  Sum  of  voltages  of  all  cycles 

%  Sum  of  thrust  for  one  cycle 
%  Kean  thrust  of  one  cycle 
%  Sum  of  mean  thrust  of  all  cycles 

■lean  frequency  of  one  cycle 
%  Sum  of  frequencies  of  all  cycles 
%  Variable  for  man  frequency 
%  Kin  free  counter 
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%  Counter  variables 

nr_cycles2=-l;         %  Counter  for  total  number  cf  cycles  in  data 
cycle_data_pts=0;     i  Number  of  data  pts  within  one  cycle 
total_data_pts2=0;     %  Number  of  data  pts  of  all  cycles 

y=-l;  %  Counter  to  determine  cycles  (25  trips  =  1  cycle) 

data_points2=length(data2 ( : , 1) ) ;  %  Counts  rows  of  data  matrix 
for  j=l :data_points2   %  Looks  at  all  data  points 
if  abs(data2 ( j , 2) )  >  5 

error  ('**  Voltage  is  outside  rangefinder  limits  (+/   5V)  **') 
end 
if  nr_cycles2  >=  0  %  Ignores  partial  data  before  1st  cycle 

cycle_volts  =  cycle_volts  +  data2  ( j-1, 2)  ;  %  Sums  cycle  volts 
cycle_thrust  =  cycle_thrust  +  data2 ( j-1, 5) ;  %  Sams  cycle  thrust 
cycle_data_pts=cycle_data_pts+l;  %  Counts  data  pts  in  cycle 
end 

I  Is  microswitch  tripped? 
if  data2 ( j , 8) >=sig_trip  4  sig_oldl<sig_trip  &  ... 

sig_old2<sig_trip  &  sig_old3<sig_trip  %  This  assumes  a  neg 
%  trip  voltage.  Change  >,  <  signs  if  tr.p  is  a  +  voltage. 
y=y+l; 
if  nr_cycles2  ==  -1 

nr_cycles2  =  0; 
end 
if  y  ==  25 

nr_cycles2=nr_cycles2+l ;   %  25  Pulses/ cycle 

y=0;  %  Resets  counter 

%  Volts 

cycle_mean_volts=cycle_volts/cycle_data_pts;  %  Cycle  mean 

data2 ( j -1 , 4 ) =cycle_mean_volts;  %  Stores  cycle  mean  volts 

volts_sum2=volts_sum2+cycle_mean_volts;  %  Sums  total  volts 

%  Thrust 

cycle_mean_thrust=cycle_thrust/cycle_data_pts;  %  Cycle  mean 

data2 ( j -1 , 6) =cycle_mean_thrust;  %  Stores  cycle  mean  thrust 

thrust_sum2=thrust_sum2+cycle_mean_thrust;  %  Sums  tot  thrust 

%  Frequency 

cycle_mean_freq=sample_rate/cycle_data_pts;  %  Cycle  mean 

data2 ( j-1 , 7) =cycle_mean_freq;  %  Stores  cycle  mean  frequency 

f req_sum2=freq_sum2+cycle_mean_freq;  %  Sums  ail  cycle  freqs 

if  data2 ( j-l,7)>0 

f req2 (x)  =  cycle_mean_f req;   %  Save  nonzero  frees 
x=x+l; 
end 

total_data_pts2=total_data_pts2+cycle_data_pts; 
%  Reset  variables 

cycle_volts=0 ;      %  Reset  cycle  voltage  for  new  cycle 
cycle_thrust=0;    %  Reset  thrust  summation  for  new  cycle 
cycle_f req=0;      %  Reset  freq  summation  for  new  cycle 
cycle_data_pts=0;  -j  Reset  nr  of  data  pts  for  new  cycle 
end 
end 

sig_oldl  =  data2(j,8);      %  Last  signal  value 
if  j>=3 

sig_old2  =  data2 ( j -1, 8) ;  %  The  data  value  before  the  last  one 
sig_old3  =  data2 ( j -2 , 8) ;  %  Three  data  points  back 
end 
end 
end 

if  k  ==  3 

data3 ( : , 1) =rawdata ( : , 2) ;  %  CH  3  is  loaded  to  datal  col  1 

data3 ( : , 8) =rawdata( : , 1) ;  %  CH  A  is  loaded  to  datal  col  8 

%  ********  Converts  data  to  voltage  and  distance  (lV=lczn)  ********** 

data3 ( :,2)=( (voltsperdiv*5 . 12) - ( (data3(:,l) )* ( vol tsperdiv* 10 .24/256) ) ) ; 

data3 ( :,9)  =  ( ( vol tsperdiv* 5 . 12 )- ( (data3 ( :,8) )* (vol tsperdiv* 10 .24/256)  )  )  ; 

data3 ( : ,3) =data3 ( : , 2) ;  %  Does  nothing 

%  --******  Converts  voltage  to  thrust  [N]  ************************** 
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value  from  controller 
value  from  controller 


(0.1  Hz  fiap 
(Thrust  runs) 


points  pno: 
trip  signal 
a  +  voltage. 


data3 ( : , 5) =  weight  *  tan (asin ( (data3 (  :  ,  3) -disp_initial ) /wire_length) ) ; 
%  ********  Computes  the  avg  dev  of  thrust  &  free  ir.  data  *********** 
%  DSO  data  variables 
if  i==100 

sig_trip  =  132; 
else  sig_trip  =  140; 
end 

sig_oldl  =  257 
sig_old2  =  257 
sig_old3  =  2  57,- 
%  Voltage  variables 
cycle_volts=0; 
cycle_mean_volts=0;   % 
volts_sum3=0;        % 
%   Thrust  variables 
cycle_thrust=0;      % 
cycle_mean_thrust=0;  % 
thrust_sum3=0 ; 
%  Frequency  variables 
cycle_mean_freq=0;    % 
freq_sum3=0; 
freq3=0; 
x=l; 

%  Counter  variables 
nr_cycles3=-l; 
cycle_data_pts=0;    % 
total_data_pts3=0; 
y=-l; 


Mir.  data 
Mir-  data 

sig_oldl ,  2 ,  and  3  are  the  3  data 
to  the  current  pt .  Use  25"  if  the 
is  neg  voltage.  Use   1  if  trip  _s 

%  Sum  of  voltage  for  one  cycle 
Mean  volts  of  one  cycle 
Sun  of  voltages  of  all  cycles 

Sum  of  thrust  for  one  cycle 
Mean  thrust  of  one  cycle 
%  Sum  of  mean  thrust  of  all  cycles 


Mean  frequency  of  one  cycle 
%  Sum  of  frequencies  of  all  cycles 
%  Variable  for  min  frequency 
%  Min  frea  counter 


Drag  run) 


%  Counter  for  total  number  of  cycles  in  data 
b  Number  of  data  pts  within  one  cycle 
:s  Number  of  data  pts  of  ail  cycles 

%  Counter  to  determine  cycles  (25  trips  =  1  cycle) 
data_points3=length(data3 ( : , 1) ) ;  %  counts  rows  of  data  matrix 
for  j=l :data_points3   %  Looks  at  all  data  points 
if  abs(data3 (j,2) )  >  5 

error  ('**  voltage  is  outside  range finder  limits  (+/-  5V)  **') 
end 
if  nr_cycles3  >=  0  %  Ignores  partial  data  before  1st  cycle 

cycle_volts  =  cycle_volts  +  data3  ( j  -1 ,  2  )  ,-  %  Sums  cycle  volts 
cycle_thrust  =  cycle_thrust  +  data3 ( j -1 , 5) ;  %  Sums  cycle  thrust 
cycle_data_pts=cycle_data_pts+l ;  %  Counts  data  pts  in  cycle 
end 


i   Is  nu.cros-.vi  ten  tripped? 
. f  data3 ( j , 8)>=sig_trip  & 


sig_oldl<sig_trip  & 


sig_old2<sig_trip  &  sig_old3<sig_trip  %  This 


signs 


f  t.i 


id  is  a 


assumes  a  neg 
-  voltage. 


^5 


Pulses /cycle 
%  Resets  coun 


;er 


%  trip  voltage.  Change  >,  < 

y=y+l; 

if  nr_cycles3  ==  -1 
nr_cycles3  =  0; 

end 

if  y  ==  25 

nr_cycles3=nr_cycles3+l; 

y=0; 

%  Volts 

cycle_mean_volts=cycle_volts/cycle_data_pts;  %   Cycle  mean 

data3 ( j-1, 4) =cycle_mean_volts;  %  Stores  cycle  mean  volts 

volts_sum3=volts_sum3+cycle_mean_volts;  %  Sums  total  volts 

%  Thrust 

cycle_mean_thrust=cycle_thrust/cycle_data_pts;    %  Cycle  mean 

data3 ( j-1 , 6) =cycle_mean_thrust;    %    Stores   cycle  mean    thrust 

thrust_sum3  =  thrust_suiri3+cycle_rnean_thrust;    %    Sums   tot    thrust 

%   Frequency 

cycle_mean_freq=sample_rate/cycle_data_pts;  %  Cycle  mean 

data3 ( j-1 , 7) =cycle_mean_freq;  %  Stores  cycle  mean  frequency 

f req_sum3=f req_sum3+cycle_mean_f req;  %  Sums  all  cycle  freqs 

if  data3 (j-l,7)>0 

f  req3  (x)    =   cycle_mean_freq;      %    Save   nonzero    :' 
x=x+l; 
end 
total_data_pts3=total_data_pts3+cycle_data_pts; 


:ecs 


sset 


.'ariables 
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cycle_volts=0;      %  Reset  cycle  voltage  for  new  cycle 
cycle_thrust=0;    %  Reset  thrust  summation  for  new  cycle 
cycle_f req=0 ;      %  Reset  free  summation  for  new  cycle 
cycle_data_pts=0;  %  Reset  nr  of  data  pts  for  new  cycle 
end 
end 

sig_oldl  =  data3 (j , 8) ;      %  Last  signal  value 
if  j>=3 

sig_old2  =  data3 ( j -1 , 8) ;  %  The  data  value  before  the  last  one 
sig_old3  =  data3 ( j -2 , 8) ;  %  Three  data  points  back 
end 
end 
end 
end 

h   Voltage  and  frequency  calculations: 

data= [ datal ; data2 ; data3 ] ,- 

total_data_pts=total_data_ptsl+total_data_pts2+total_data_pts3 ; 

vol ts_sum=vol ts_suml+vol t s_sum2  +vol ts_sum3  ; 

f req_sum=  f req_suml  +  f req_sum2+ f  req_sum3 ; 

thrust_sum=thrust_suml+thrust_sum2+thrust_sum3 ,- 

nr_cycles=nr_cyclesl+nr_cycles2+nr_cycles3; 

freq=[freql  freq2  freq3]; 

clear  datal 

clear  data2 

clear  data3 

volts_mean=volts_sum/nr_cycles; 

if  volts_mean==0  %  No  or  very  slow  flap  case,  still  need  data 

volts_mean=mean (data ( : , 3 ) ) ; 
end 

inner_suml=0;   %  Inner  summation  for  adev  of  volts 
inner_sum2=0 ;   %  Inner  summation  for  adev  of  freq 
freq_mean=f req_sum/nr_cycles;  %  Mean  freq  of  cycles 
max_freq=max ( f req) ;  \   Max  cycle  freq  of  all  cycles 
min_freq=min( f req) ;  %  Min  cycle  free  of  all  cycles 
data_points=length(data( : ,1) > ; 
for  m=l :data_points 

if  data(m,6)~=0 

inner_suml=inner_suml  +  abs (data(m, 4) -volts_mean) ; 
inner_sum2  =  inner_sum2  +  abs (data (m, 7) -f req_mean) ; 

end 
end 

adev_f req=inner_sum2/nr_cycles;  %  Avg  dev  of  frequency 
%  Thrust 

thrust_mean=weight*tan(asin( (volts_mean-disp_initial) /wire_length) ) ;  %  Mean  thrust  of 
cycles 

%  **'"  There  are  4  sources  of  thrust  error: 

%  ***  1st  one  is  due  to  cyclic  nature  of  flapping  wing  thrust:  ***** 

adev_volts=inner_suml/nr_cycles;  %  adev  fm  voltage  data 

%  *  *  *  2nd  source  of  error  -  laser  rangefinder  error  [mmj  *********** 

if  abs (volt s_mean)  <  3.5 

error_laser  =  0.1  +  0.002  *  abs (10  *  volts_mean) ; 
else 

error_laser  =  0.25  +  0.005  *  absdO  *  volts_mean)  ; 
end 

error_laser  =  error_laser/10;   %  Convert  range  error  to  volts  error 
%  +**  Oscilloscope  volt  max  error=resolution  of  DSO  2102- 3rd  error  * 
error_DSO  =  (voltsperdiv*5 . 12) /128; 

%  ***  Convert  errors  fm  voltage  to  thrust  [N]  error  **************** 
adev_thrustl=weight*tan(asin( (adev_volts) /wire_length) ) ; 
error_thrust2=weight*tan(asin( (error_laser) /wire_length) ) ; 
error_thrust3=weight*tan (asin( (error_DSO) /wire_length) ) ; 

s  ***  Computes  error  due  to  v.'eight  and  wire  length  -  4th  error:  ********* 
error_thrust4  =weight_max*tan (asin( (volts_mean-disp_initial) /wire_length_min) ) - 
thrust_mean; 

error  thrust=adev_thrustl+error_thrust2+error_thrust3+error_thrust4; 
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%  *****,,***  *******  saves  ail  data  co  filename.wri  ************  ****** 
%   Used  for  trouble  shooting  the  program 

%filesave2= (  [  'H:  \cocs\Resuits'-  '  , num2str  (i )  ,  filename,  '  .v;ri  '  ] )  ; 
%  i ic2  =  f open ( f i lesave2 ,  '  w'  !  ; 
%fprintf (fid2,  [numzstr ( i ; , filename! j ; 
%fpnntf  !f  ic2,  '  \rCH  B\tCH  B\tCH  B\tCH  B\tCH  B\t' )  ; 
Sfprintf (fid2, 'CH  3\tCH  3\tCH  A\tCH  A\n'!; 

%fprir.tf  !fid2,  'Data\tVoltage\ tAfosol  Volt\tCycle  Mean\tThrust\t '  )  ,- 
%fprintf  if  id2,  'Cycle  Mean\t?req  Kean\tData\-Vcltage\n' i  ; 
%f  print  f  !f  ic2,  '\t  [V]  ,  [cm]  \t  [Vj ,  [cm]  \t  [V]  \t  [N]  \C  [N]  \t ' )  ; 
Sfprintf if id2, ' [Hz] \t\t [V]  \n' ) ; 

%fprintf ! f id2, ' \n\tMean  thrust  =\t\t\c  %10.3f ', thrust jnean) ; 
%fprintf  (fid2,  '  \r.\tThrust  Ave  Dev  Error  =\t\t\t  %10. 8f , adev_thrustl) ; 
%fprintf  if  id2  ,  '  \n\tLaser  Error  =\t\t\t  %10.8f  ,error_thrust2)  ; 
"ofprmtf  (f  id2,  '  \n\tOscilioscope  Error  =\t\t\t  %10.Sf  ,error_thrust3)  ; 
%fprintf (fid2, '\n\tWt  Wire  Error  =\t\t\t  %10  .  Sf  '  ,  error_thrust4}  ,- 
%f print f (fid2, ' \n\tTotal  Error  =\t\t\t  %10 . Sf ' , error_thrust i ; 
%fprintf (fid2, ' \n\t\t\t\tMean  frequency  =\t\t  %4.2f\n' 
%    freq_mean) ; 

%fprintf  (fid2, '\t\t\t\tAvg  Dev  =\tv-t  %4.2f  \n'  ,adev_£req) ; 
%fprintf (fid2, ' \t\t\t\tMin  Freq  =\t\t  %4 .2f \n' ,min_freq) ; 
%fprintf (f id2 , ' \t\t\t\tMax  Freq  =\t\t  %4 . 2f \n\n' ,rrax_f req) ; 

%fprintf (fid2, '%3.0f\t  %4.3f\t  %5.4f\t  %5.4f\t  %S.5f\t  %8.5f\t  %5.2f\t  %3.0f\t 
%4 .3f\n' ,data. ' ) ; 
%fclose(f id2  i ; 

%  »********+*,+  saves  summary  of  data  to  Mean_Thrust.wri  *********** 
fid3  =  fopen( 'H:\docs\Results\Mean_Thrust.wri' , 'a' ) ; 
fprintf (fid3, [num2str (i) , filename] ) ; 

fprintf (fid3, '\t\t\t%10.8f \t\t%10.8f\t%10.8f \t%10.8f\t%10.8f \t%10.8f\t%4.2f\t%4.2f\t%4.2f 
\t%4.2f\n' , . . . 

thrust_mean, error_thrust , adev_thrustl, error_thrust2 , error_thrust3 , . . . 

error_thrust4 , freq_mean, adev_f req,min_f req,max_f req) ; 
fclose(fid3) ; 
end 
disp( ['*********  Program  Complete  ***********']) 
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APPENDIX  B.  THREE  DIMENSIONAL  PLOTS  OF  EXPERIMENTAL  DATA 
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Figure  36.  3-D  Plot  of  Low  Stiffness  Model's  Thrust  Performance. 
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Figure  37.  Contour  Plot  of  Low  Stiffness  Model's  Thrust  Performance. 
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Medium  Stiff  Mode! 
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Figure  38.  3-D  Plot  of  Medium  Stiffness  Model's  Thrust  Performance. 
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Figure  39.  Contour  Plot  of  Medium  Stiffness  Model's  Thrust  Performance. 
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Nigh  Stiff  Model 
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Figure  40.  3-D  Plot  of  High  Stiffness  Model's  Thrust  Performance. 
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Figure  41.  Contour  Plot  of  Medium  Stiffness  Model's  Thrust  Performance. 
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Articulated  Pitch  Mods! 
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Figure  42.  3-D  Plot  of  Articulated  Pitch  Model's  Thrust  Performance. 
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Figure  43.  Contour  Plot  of  Medium  Stiffness  Model's  Thrust  Performance. 
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